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Abstract 



A measurement of the tt production cross section in proton-proton collisions at a 
centre-of-mass energy of 7 TeV has been performed at the LHC with the CMS de- 
tector. The analysis uses a data sample corresponding to an integrated luminosity of 
36 pb _1 and is based on the reconstruction of the final state with one isolated, high 
transverse-momentum electron or muon and three or more hadronic jets. The kine- 
matic properties of the events are used to separate the tt signal from W+jets and QCD 
multijet background events. The measured cross section is 173^2 (stat. + syst.) pb, 
consistent with standard model expectations. 
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1 Introduction 

The top quark occupies a unique position within the standard model. With a mass roughly 
that of a tungsten atom, it is the only quark heavy enough to decay before forming bound 
states with other quarks. Its large mass has inspired numerous theoretical models in which the 
top quark plays a special role in the generation of mass or in the physics of new, undiscovered 
particles. The top quark often acts as either a direct contributor to new physics or an important 
background in new-particle searches in these models. 

In hadron colliders, top-quark production is dominated by the production of tt pairs Q. At 
the Tevatron, where the top quark was discovered in 1995 [2, 3], tt pairs are predominantly 
produced through quark-antiquark annihilation. In contrast, the tt production mechanism at 
the Large Hadron Collider (LHC) [4] is expected to be dominated by the gluon fusion process. 
Measurements of the tt cross section at the LHC can provide important tests of our understand- 
ing of the top-quark production mechanism and can also be used in searches for new physics. 

In the standard model, a top quark decays nearly 100% of the time to a W boson and a b quark. 
The decay of a tt pair is categorized by the decay of the W bosons produced by the pair. Thus 
the channel in which both W bosons decay to leptons is referred to as the "dilepton" chan- 
nel, and the channel in which one W decays to leptons and the other to quark jets is the 
"lepton+jets" channel. The channel in which both W bosons decay to jets is called the "all 
hadronic" channel. A further categorization of the decays is made by specifying the flavour 
of the charged lepton(s) produced from the W decays. For the purposes of this paper, the 
"lepton+jets" channel refers only to decays in which the charged lepton is either an electron or 
a muon. 

The next-to-leading-order (NLO) corrections to the top-quark pair production cross section at 
hadron colliders have been calculated both with the full top-quark spin dependence [5 . 6] and 
without this dependence 0151 . A complete analytic result for the NLO partonic cross section 
has only been published recently [9[. Approximations of a full next-to-next-to-leading-order 
(NNLO) calculation have also been obtained by various groups IIT0HT31 and references therein]. 
The tt production cross section at the LHC has been previously measured in the dilepton chan- 
nel by the Compact Muon Solenoid (CMS) experiment [14) and in the combined dilepton and 
lepton+jets channels by the ATLAS experiment [15J. These measurements agree with recent 
NLO and with approximate NNLO calculations. 

In this paper we present a measurement of the cross section for tt production in proton-proton 
collisions at y/s = 7 TeV with the CMS detector, using the electron +jets and muon+jets final 
states. Although there are two jets from hadronization of the b quarks in these final states, in 
this analysis no requirement is made on the presence of b jets. Instead, the kinematic properties 
of the events are used to select the tf signal. It is important to measure the tf cross section both 
with and without a requirement on the presence of b jets, as new physics could contribute 
differently in each case. 

A brief overview of the CMS detector is provided in Section [2] of this paper, followed by a 
discussion of event reconstruction procedures in Section |3j The selection criteria applied to 
the data are described in Section |4] and the processes used to simulate signal and background 
events are described in SectionlU Section|6]details the method used to extract a measurement of 
the tt cross section from the selected events, as well as the calculation of statistical and system- 
atic uncertainties on the result. Section|7]summarizes the result and compares the measurement 
with recent perturbative QCD calculations. 
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3 Event Reconstruction 



2 CMS Detector 

The central feature of the CMS apparatus is a 3.8 T magnetic field produced by a superconduct- 
ing solenoid of 6 m internal diameter. Within the field volume are the silicon pixel and strip 
trackers, the crystal electromagnetic calorimeter (ECAL), and the hadron calorimeter (HCAL). 
The inner tracker measures charged particles within the range \rj\ < 2.5, where t] indicates 
detector pseudorapidity. It consists of 1440 silicon pixel and 15 148 silicon strip detector mod- 
ules and is located within the axial magnetic field. It provides an impact parameter resolution 
of ~ 15 fim and a transverse momentum (pj) resolution of ~ 1.5% for 100 GeV particles. The 
ECAL consists of nearly 76 000 lead tungstate crystals that provide coverage in pseudorapid- 
ity of 1 77 1 < 1.48 in the ECAL barrel region and 1.48 < \t}\ < 3.0 in the two endcaps. A 
preshower detector consisting of two planes of silicon sensors interleaved with a total of 3 ra- 
diation lengths of lead is located in front of the endcaps. The ECAL energy resolution is 3% or 
better for the range of electron energies relevant for this analysis. 

The HCAL is composed of layers of plastic scintillator within a brass/ stainless steel absorber, 
covering the region \tj\ < 3.0. A calorimeter composed of quartz fibres embedded in a steel 
absorber extends the forward HCAL coverage beyond the solenoid volume, to \t] \ < 5.0. In the 
region \t]\ < 1.74, the HCAL cells have widths of 0.087 in both pseudorapidity and azimuth 
(</>). In the (t],(p) plane, for \rj\ < 1.48, the HCAL cells match the underlying 5x5 ECAL crystal 
arrays to form calorimeter towers projecting radially outwards from close to the nominal inter- 
action point. At larger values of \rj\, the size of the towers increases, and the matching ECAL 
arrays contain fewer crystals. 

Muons are measured in the pseudorapidity window \tj\ < 2.4, with detection planes made of 
drift tubes, cathode strip chambers, and resistive plate chambers. Matching the muons to the 
tracks measured in the silicon tracker results in a transverse-momentum resolution between 
1% and 5% for pj values up to 1 TeV. A two-tier trigger system selects the most interesting pp 
collision events for use in physics analyses. A more detailed description of the CMS detector 
can be found elsewhere [16|. 

3 Event Reconstruction 

The reconstruction of electrons [flTl [T8| uses information from the pixel detector, the silicon 
strip tracker, and the electromagnetic calorimeter. The amount of material before the ECAL in 
the CMS detector ranges from 0.4 to 1.6 radiation lengths, depending on t], and an electron may 
lose a considerable fraction of its energy through bremsstrahlung in passing through this mate- 
rial. The energy deposited in the calorimeter may be spread over a wide range in <p compared 
to the initial direction of the electron. To account for this in the reconstruction of electron can- 
didates, clusters of calorimeter energy deposition ("superclusters") from a narrow fixed range 
in 77 and a variable range in (p are formed. Starting from these superclusters, corresponding 
hits in at least two layers of the pixel tracker capable of acting as seeds for electron trajectory 
candidates are required. Energy loss through bremsstrahlung leads to non-Gaussian contribu- 
tions to fluctuations in the calorimeter and tracking measurements. Therefore, the seeding and 
building of tracks is done using dedicated algorithms designed to handle these fluctuations. 
The final fit of the trajectories relies on a Gaussian sum filter that is a non-linear generalization 
of the Kalman filter with weighted sums of Gaussians instead of a single Gaussian. 

Multiple reconstruction algorithms exist to identify muon candidates in CMS from hits in the 
silicon tracking system and signals in the muon system |ll9ll20| . Since muons are typically the 
only particles reaching the muon chambers, the "stand-alone muon" reconstruction algorithm 
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starts from track segments detected in the innermost layers of the muon chambers. Additional 
hits in surrounding layers are then added. Finally muon tracks are propagated back to the in- 
teraction point. A second reconstruction algorithm ("tracker muons") starts from tracks found 
in the tracking system and then associates them with compatible signals in the calorimeters and 
the muon chambers. A third algorithm ("global muons") starts from stand-alone muons and 
then makes a global fit to the consistent hits in the tracker and the muon system. Muons used 
in this analysis are required to be identified by both the tracker muon and the global muon 
algorithms. 

Hadronic jets are reconstructed from individual particles whose identities and energies have 
been determined via a particle-flow technique pl| that combines information from all subde- 
tectors: charged tracks in the tracker and energy deposits in the electromagnetic and hadronic 
calorimeters, as well as signals in the preshower detector and the muon system. The energy 
calibration is performed separately for each particle type. All particles found by the particle- 
flow algorithm are clustered into particle-flow jets f22ll231 by using each particle's direction at 
the interaction vertex and the anti-/cx jet clustering algorithm [24J with the distance parameter 
R = 0.5, as implemented in FASTjET version 2.4 1125 1 . 

Since most jet constituents are identified and reconstructed with nearly the correct energy by 
the particle-flow algorithm, only small residual jet energy corrections must be applied to each 
jet. These corrections are between 5% and 10% of the jet energy and were obtained as a function 
of p T and ?/ from the GEANT4-based CMS Monte Carlo simulation (v. 9.2 RevOl) |26[ and early 
collision data. The factors also include corrections for small discrepancies observed between 
the simulation and the data. 

The missing transverse energy (Ej) is defined as the magnitude of the negative vector sum of 
the transverse energies (Ej) of all the particles found by the particle-flow algorithm. A decay 
of a tt pair via the lepton+jets channel is expected to exhibit significant missing tranverse en- 
ergy because of the undetected neutrino from the leptonically decaying W. Distributions of this 
variable are used in likelihood fits to measure the tt signal and to distinguish it from various 
backgrounds, as discussed in the following sections. 

4 Data Set and Event Selection 

The data discussed in this paper were collected in the period April to November 2010 from 
proton-proton collisions at y/s = 7 TeV and correspond to an integrated luminosity of 
36 ± 1 pb -1 ||27ll28|. The trigger providing the data sample used in this analysis is based on 
the presence of at least one charged lepton, either an electron or a muon. Because the peak in- 
stantaneous luminosity increased throughout the data-taking period, the minimum transverse 
momentum pj of the muon required in the trigger ranged from 9 to 15 GeV, and the minimum 
Ej required in the trigger for electrons similarly ranged from 10 to 22 GeV. This data sample is 
used for the selection of the signal region and the selection of signal-depleted control regions 
used for studies related to background processes. Trigger efficiencies are determined from the 
data using Z-boson events and then corrected for the differences in the efficiencies between 
the Z-boson events and tt events, using the simulated samples described below. Events are 
required to have at least one primary pp interaction vertex, where vertices are identified by 
applying an adaptive fit to clusters of reconstructed tracks |29fl. The primary vertex must be 
within ±24 cm of the nominal interaction point in the direction along the proton beams. The 
distance between the primary vertex and the nominal interaction point in the plane perpendic- 
ular to the beam direction must be less than 2 cm. 
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4 Data Set and Event Selection 



In the event selection for the electron +jets channel, at least one electron with transverse en- 
ergy greater than 30 GeV and \t]\ less than 2.5 is required. Electrons from the transition region 
between the ECAL barrel and endcap, 1.44 < \rj sc \ < 1.57, are excluded, where t] sc is the pseu- 
dorapidity of the ECAL supercluster. The energy of the HCAL cell mapped to the supercluster 
must be less than 2.5% of the total calorimeter energy associated with the supercluster. Ad- 
ditional requirements are made on the shower shape and the angular separation between the 
ECAL supercluster and the matching track. Electron tracks must extrapolate to within 0.02 cm 
of the interaction vertex in the plane perpendicular to the proton beams and to within 1 cm 
in the direction along the beams. Electron candidates that lack hits in the inner layers of the 
tracking system are assumed to be the product of photon conversions and are discarded. 

Since the electron from the W boson in a top-quark decay is expected to be isolated from other 
high-pj particles in the event, electrons are required to have a relative isolation (7 re i) smaller 
than 0.1, where relative isolation is defined as the scalar sum of the transverse momenta of 
tracks with pj > 1 GeV and all calorimeter energy in a cone of AR = y/ (A</>) 2 + (At]) 2 < 0.3 
around the electron, divided by the electron pj. Here Acp (At]) is the difference in azimuthal 
angle (pseudorapidity) between the electron and the track or calorimeter cell. Contributions to 
the sum due to the electron itself are removed. Events containing multiple electron candidates 
are rejected if any combined dielectron invariant mass lies within 15 GeV of the Z-boson mass. 

In the event selection for the muon+jets channel, muons are required to have at least a mini- 
mum number of hits in both the silicon tracking system and the muon chambers. The muon 
must have transverse momentum greater than 20 GeV and must lie within the muon trigger 
acceptance < 2.1). The muon must be separated from any selected jet (defined below) in 
the event by AR > 0.3. Muons must extrapolate to within 0.02 cm of the interaction vertex 
in the plane perpendicular to the beams and to within 1 cm in the direction along the beams. 
The muon from the W boson in a top decay is also expected to be isolated from other high-px 
particles in the event, and thus muons are required to have relative isolation smaller than 0.05, 
where the muon relative isolation is defined analogously to the electron J rel . Exactly one muon 
passing all these criteria is required. Events containing a separate, more loosely-defined muon, 
as well as the highly-energetic isolated muon described above, are rejected. 

Z-boson events are used to study lepton trigger, identification, and isolation efficiencies in 
the data. A high purity Z-boson sample is extracted from data by requiring two oppositely- 
charged like-flavour leptons with a combined invariant mass within 15 GeV of the Z-boson 
mass. The two-electron final state suffers from background contamination due to hadronic 
jets misidentified as electrons. This contamination is modeled using events containing a pair 
of like-charged electrons whose invariant mass falls within the Z-boson mass window. The 
identification efficiency for isolated electrons is determined to be 0.75, with combined statistical 
and systematic uncertainty (stat.+syst.) of ±0.01. The trigger efficiency for such electrons is 
measured to be 0.982 ± 0.001 (stat). Using Z — > events, the efficiency of finding an isolated 
muon is measured to be 0.880 ± 0.002 (stat. + syst), and the efficiency for triggering on muons 
passing all selection cuts is found to be 0.922 ± 0.002 (stat.). 

Selected jets are required to have a jet-energy-scale-corrected px > 30 GeV, \tj\ < 2.4, and must 
be separated by AR > 0.3 from isolated electrons and AR > 0.1 from isolated muons in order 
to remove jet candidates produced by charged leptons. Both the muon+jets and electron +jets 
analyses require that events contain at least three jets. Selected events are grouped into sub- 
samples based on their jet multiplicity, so that events containing exactly three jets are separated 
from those containing four or more jets. A requirement on is not included in the event se- 
lection, as fits to the distribution of this observable are used to separate the tt signal from the 
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backgrounds. 

Dilepton tf decays are removed by discarding events that contain both a high-px electron and 
a high-px muon. In the electron+jets event selection, events containing any muon with p x > 
10 GeV, \rj\ < 2.5, and relative isolation < 0.2 are rejected. In the muon+jets event selection, 
events must contain no electron candidates with Ej > 15 GeV and 7 re i < 0.2. 

5 Simulation of Signal and Background Events 

We determine the efficiency for selecting lepton+jets signal events by using a simulated tf 
event sample. We perform the simulation of signal tt events using MadGraph (v. 4.4.12) [30 1, 
where the top-antitop pairs are generated accompanied by up to three additional hard jets in 
the matrix-element calculation. The factorization and renormalization scales Q 2 for tt are both 
set to 

Q 2 = m 2 op + Zp 2 , (1) 

where m t0 p is the top-quark mass and Ep T is the sum of the squared transverse momenta of all 
accompanying hard jets in the event. MadGraph is also used to generate background events 
from electroweak production of single top quarks, the production of leptonically-decaying W 
and Z bosons in association with up to four extra jets, and photon+jets processes. The factor- 
ization and renormalization scales in these events are set in the same manner as for tt events, 
with the appropriate boson mass replacing the top-quark mass. The parton-level results gen- 
erated by MadGraph are next processed with Pythia (v. 6.420) |3"T1 to provide showering 
of the generated particles. Shower matching is done by applying the MLM prescription 11321 . 
Events are then passed through the GEANT4-based CMS detector simulation. 

In addition to the simulated events generated with MadGraph, several QCD multijet samples 
have been produced using Pythia. Samples enriched in electrons (muons) are used to provide 
a preliminary estimate on the QCD background in the electron+jets (muon+jets) channel. The 
final QCD background contribution is taken directly from data, as is explained in Section |6j 

The tf NLO production cross section has been calculated as a t i = 157+24 pb, using MCFM 11331 
[34j. For both tf and single-top-quark production (described below), renormalization and fac- 
torization scales were set to Q 2 = (172.5 GeV) 2 . The uncertainty in the cross section due to 
uncertainties in these scales is determined by varying the scales by factors of 4 and 0.25 around 
their nominal values. Contributions to the cross section uncertainty from the parton distribu- 
tion functions (PDF) and the value of as are determined following the procedures from the 
MSTW2008 B35|, CTEQ6.6 1351, and NNPDF2.0 |37| sets. The uncertainties are then combined 
according to the PDF4LHC prescriptions II38II . 

The t-channel single-top-quark NLO cross section (multiplied by the leptonic branching frac- 
tion of the W boson) has been determined as a t = 21.0+}; J pb using MCFM I551I59H4H, where 
the uncertainty is defined similarly as for tf production. The inclusive single-top-quark asso- 
ciated production (tW) NLO cross section of (7 tw = 10.6 ± 0.8 pb Il40l has been used. Both 
cross sections include the production of single top and single antitop quarks. The s-channel 
single-top production cross section is small compared to the t-channel and tW production cross 
sections and is treated as negligible in this analysis. 

The NNLO production cross section for W bosons decaying into leptons has been determined 
to be c%_^„ = 31.3 ± 1.6 nb using FEWZ |42| . Its uncertainty was determined in a similar 
manner as for top-quark pair production. Finally, the Drell-Yan dilepton (££) production cross 
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6 Cross-Section Measurement 



section at NNLO has been calculated using FEWZ as cr z/r ^ u (m u > 20 GeV) = 5.00 ± 0.27 nb 
and oz/y*-^u( m U > 50 GeV) = 3.05 ± 0.13 nb. Backgrounds due to diboson production have 
been ignored given their relatively small expected contribution to the lepton+jets event yield. 

In a simulated sample of tt events in which all top-quark decay modes are included, the 
electron+jets selection efficiency is found to be 5.7%, while the muon+jets selection efficiency 
is 7.2%. The selected simulated signal events in each mode are dominated by tt decays to 
electron+jets and muon+jets, respectively, although tt decays containing tau leptons also con- 
tribute. Table[l]and Fig.[I]give the observed numbers of events in both channels after applying 
the event selection procedures described above to the pp collision data set. The numbers of 
events predicted by the simulation for the different physics processes are also listed. Predicted 
yields are calculated by multiplying selection efficiencies for each process, as determined from 
simulation, by the appropriate NLO or NNLO cross sections and the total integrated lumi- 
nosity of 36 pb _1 . As shown in Fig. [lj the fractional contribution to the event yield from tt 
signal events is negligible in events with zero jets, but dominates as jet multiplicity increases. 
More events are observed in the data than predicted by the simulation, indicating that either 
the signal cross section or the background cross sections, or possibly both, are larger than ex- 
pected. Our method for determining the number of signal and background events in the data 
is described in the next section. 

Table 1: Predicted event yields for the electron+jets and muon+jets event selections. The event 
yields from the simulation are normalized to an integrated luminosity of 36 pb -1 . The quoted 
uncertainties account for the limited number of simulated events, the uncertainty on the cal- 
culated cross section (where available), the uncertainties on the lepton selection and trigger 
efficiency correction factors, and a 4% uncertainty on the luminosity. The penultimate row lists 
the totals from the simulation, and the last row shows the number of observed events. 





electron+jets 

Njets = 3 N jets > 4 


muon+jets 

Njet = 3 Njets > 4 


tt 

single top 
W+jets 
Z+jets 
QCD 


157 ± 25 168 ± 27 
22 + 1 9 + 1 

374 ±27 94 + 7 
66 + 5 15 ± 1 

314 ± 19 53 ± 8 


197 ± 31 211 ± 33 
30 ± 1 11 ± 1 

486 ± 34 115 ± 9 
46 ± 3 11 ± 1 
49 ± 3 9 ± 1 


sum (simulated events) 


934 ± 55 339 ± 32 


807 ± 53 358 ± 37 


observed in data 


1183 428 


1064 423 



6 Cross-Section Measurement 
6.1 Analysis Method 

We measure the tt cross section in a data sample consisting of events with leptons and jets in 
the final state. Since the event yields for the background processes can be difficult to estimate 
purely from simulations, the kinematic properties of the events are used to separate signal from 
background. It would be natural to require four or more jets for the selection of tt events in the 
lepton+jets channel, owing to the four final-state quarks present in these decays and because 
the number of background events from W/Z+jets and QCD multijet events decreases with 
increasing jet multiplicity. However, it is also useful to include events with only three jets in 
the selection. In addition to improving the overall signal efficiency, the inclusion of three-jet 
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Figure 1: Observed numbers of events from data and simulation as a function of jet multiplicity 
in the (left) electron+jets and (right) muon+jets selected samples. Yields are calculated after 
applying the respective event selections and omitting any requirement on the number of jets. 

events constrains the QCD and W+jet background normalization when a simultaneous fit to 
the three and four-or-more jet samples is performed. 

The predicted jet-multiplicity distribution and the predicted ratio between events with three 
jets and events with four or more jets for the different processes are used to simultaneously fit 
the fraction of tt events and the contamination from background processes. Kinematic variables 
whose shapes are different for the different processes are used to separate the backgrounds 
from the signal. After a number of variables and combinations were tested, the variable M3 was 
chosen to separate tt events from background events in the four-or-more-jet sample. This vari- 
able is defined as the invariant mass of the combination of the three jets with the largest vecto- 
rially summed transverse momentum. It approximates the mass of the hadronically-decaying 
top quark and thus provides good separation power. The three-jet sample is dominated by 
W+jets events and QCD multijet events. For the three-jet sample, a variable that is well suited 
for the discrimination of QCD multijet events from the other processes is needed. In contrast 
to processes with W bosons, QCD processes exhibit only small amounts of missing transverse 
energy, mostly because of mismeasured jets rather than the presence of neutrinos. Therefore 
#t was chosen as the discriminating variable to separate QCD events from W+jets and tt signal 
events in the three-jet sample. 

The #t and M3 distributions from the observed data sample are fit simultaneously to obtain 
the contributions of the signal and the main background processes. We use a binned likelihood 
fit, where the number of expected events jij [i] in each bin i of the distribution for the variable 
of choice j (either or M3) is compared to the observed number of events in this bin. The 
number of expected events in bin i is given by: 



( 2 ) 

k 

where is the binned contribution (called "template" in the following) for variable j and 

process k. The fit parameters fa are the ratio of the measured (c^) and predicted (of* ) cross 
sections for process k: 

h = 0) 

Here, k denotes all processes that are taken into account, namely tt, W+jets, Z+jets, single- 
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6 Cross-Section Measurement 



top-quark decays, and QCD. Since negative fa values are unphysical, we do not allow fa to 
become smaller than zero in the fit. The templates a» are normalized to the corresponding 
predictions for process k for events with three jets and events with four or more jets. Because 
of the normalization of the fit templates to the prediction, the fitted value of fa can be directly 
interpreted as the scale factor one has to apply to the predicted cross section of a given process 
k to derive the measured cross section. 

Templates for the tt, single-top, W+jets, and Z+jets processes are derived from the simulation, 
while the QCD multijet template is derived from data, using a method that will be described 
later. The shapes of the M3 and fyj distributions of single-top-quark events are very similar to 
the distributions of tt events. Because of this similarity and because of the very small number 
of expected single-top-quark events, an unconstrained fit of the single-top-quark contribution 
is not possible. However, since the single-top-quark production process is theoretically well 
understood, the number of such events can be estimated from simulations. Therefore the fit 
parameter for single top is not left to float freely, but is instead subject to a Gaussian constraint 
with a mean of 1.0 and width of 0.3. The uncertainty on this value is assigned according to the 
expected precision of initial single-top cross-section measurements in CMS [43 1. In addition, 
the ratio between the W+jets and Z+jets cross sections is constrained to be within 30% of the 
expectation from theory, where the constraint width is set by the uncertainty in the NLO cross 
sections [44J. These constraints are inserted by multiplying the likelihood function used in the 
fit by Gaussian terms of mean value 1.0 and widths corresponding to the uncertainties on the 
respective constraints. The same fa parameters are used for both jet-multiplicity bins. 

A Neyman construction [45J with central intervals and a maximum-likelihood estimate of the 
tt cross section as test statistic was chosen to obtain the confidence interval for the tt cross 
section. For this purpose pseudo-experiments are performed in which the number of events 
from the different processes are chosen randomly around the values predicted by simulations 
within appropriate uncertainties. Specifically this is done by randomly choosing, for each back- 
ground process, an input value for fa from a normal distribution with mean value 1.0 and a 
width of 30% for W+jets, Z+jets, and single top. Since the properties of QCD multijet events 
are more difficult to calculate, a more conservative uncertainty of 50% is used for this back- 
ground. The templates for the different processes are then scaled with the corresponding fa 
values and summed together, generating a pseudo-data distribution for in the three-jet bin 
and a pseudo-data distribution for M3 in the inclusive four-jet bin. To simulate statistical fluc- 
tuations we then vary the contents of each bin of the two distributions using Poisson statistics. 
A maximum-likelihood fit to the templates is then performed on the distributions. This proce- 
dure yields one signal fit result jS|* for each pseudo-experiment and provides a measurement 
of any possible bias between the input and fitted values. We vary the input value between 
0.0 and 3.0 in steps of 0.2, and perform 50 000 pseudo-experiments for each value. Each set of 
pseudo-experiments gives a distribution of the fitted values j6{j . For each input value we deter- 
mine the median and the 68% and 95% quantiles of the corresponding jS^ distribution and use 
these values for the estimation of the central values and for the construction of confidence belts, 
respectively. From this confidence-belt construction the tt cross section result corresponding to 
the jf?!* measured in data can be extracted together with its total uncertainty.By construction, 
this treatment the correct coverage probability. 

6.2 Systematic Uncertainties 

In general, the presence of a systematic uncertainty affects both the number of selected events 
and the shape of the investigated discriminating observables, resulting in modified distribu- 
tions a. s J st for the different processes. In order to estimate the effect of different sources of 
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systematic uncertainties we construct modified templates and draw the pseudo-data used for 
the statistical calculation from them. For each source of systematic uncertainty u, two templates 
a"^ +1 and oH^T , corresponding to variations of ±1 standard deviation (±lc) of the specific sys- 
tematic uncertainty are used. Both templates are normalized to the altered event yields for each 
specific systematic uncertainty u, and thus account for both changes in event rates and changes 
in parameter distributions. These ±lc templates are derived either by altering the nominal 
samples as described in the following sections or from dedicated simulations. The modified a" k 
used for drawing the pseudo-data can then be constructed from the nominal template ota and 
the ol 1 -^ 1 templates. Therefore for each uncertainty u a strength parameter 5 U is introduced, and 
a^ st is defined as a linear interpolation: 

4 ySt H = cc jk \i] +Y:\5 U \ ■ (4' si ^[i] - cc jk [t\) . (4) 

u 

Here, u runs over all sources of systematic uncertainties and ocj^ 1 [i] is the prediction for bin 
i of distribution j of process k affected by +lcr or —la of uncertainty u. Random numbers 
following a Gaussian distribution with a mean of zero and unit width are used as values for 
each S u , with the strengths determined separately for each pseudo-experiment. The nominal 
template is reproduced for 3 U = 0, while the two altered templates correspond to S u = +1 and 
S u = —1. For all other values of 5 U , the desired mixture of the nominal and shifted templates is 
obtained. 



ation is cut 
has a value 



In order to prevent unphysical negative predictions for a process, the linear interpo 
off at a bin content of zero, i.e., whenever a bin of oc^ st calculated according to Eq. |i 

less than zero, zero is used instead. The definition of jSS* remains unchanged, meaning that the 
original templates without uncertainties are employed for fitting the pseudo-data distributions. 

The influence on our measurement due to the imperfect knowledge of the jet energy scale (JES) 
is estimated by simultaneously varying all jet four-momenta either by +\a or by —lcr of the 
energy scale uncertainties, which are functions of the jet t] and pj. These uncertainties on the 
particle-flow jet energy scale are typically about 3%, as shown in Ref. 11461 . In addition, a con- 
stant 1.5% uncertainty due to changes in calorimeter calibrations and a px-dependent uncer- 
tainty of 1.32 GeV/ pj due to multiple collisions in the same event ("pileup") are applied. For 
jets that can be matched to a b quark at the parton level, we assign an additional uncertainty of 
2% if the jet lies within \ tj\ < 2.0 and has pj between 50 GeV and 200 GeV, and a 3% uncertainty 
otherwise. This uncertainty accounts for observed response differences for b jets generated in 
Pythia and those in Herwig |47|. The overall uncertainty is determined by adding all of the 
individual uncertainties in quadrature. 

Jet asymmetry measurements suggest that the jet energy resolutions (JER) in data are about 
10% worse than in the simulation [48 1. The uncertainty on this measurement is also about 10%. 
To account for this, jets in the simulated samples are altered so that their resolutions match 
those measured in data, and the effect is propagated to the calculation of #t- The impact of this 
uncertainty on our measurement is determined by evaluating the change in cross section when 
simulated jet resolutions are widened by 0% or 20%, rather than the default 10%. 

The corrections in jet energy scale and resolution described above are used to vary the missing 
transverse energy according to variations in clustered jet energy. In order to also account for 
the uncertainty of unclustered energy in Uj, the amount of unclustered energy contributing to 
#x is shifted by ±10%. However, the impact of the variation of the unclustered energy on the 
measurement is found to be negligible. 
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6 Cross-Section Measurement 



Adjusted simulated samples are used to evaluate the dominant systematic uncertainites in the 
cross section measurements. Two simulated samples of tt events are available to estimate the 
systematic uncertainty induced by the lack of accurate knowledge of the amount of QCD initial- 
state (ISR) and final-state (FSR) radiation. For one of these samples, the amount of ISR and FSR 
has been increased, while less ISR and FSR is assumed in the other sample. 

The impact of uncertainties in the factorization scales on the cross section measurement is es- 
timated by varying the scales in each of the samples by factors of 0.25 and 4.0 with respect 
to their default values. The W+jets and Z+jets processes are treated as being correlated, and 
their respective factorization scales are shifted either down or up simultaneously, while the tt 
scale is considered to be uncorrelated and is shifted independently. The impact of a variation 
of the shower matching threshold is investigated by varying the matching thresholds for the 
three processes by factors of 0.5 and 2.0 compared to the default thresholds. Again, W+jets 
and Z+jets processes are treated as fully correlated and are varied simultaneously, while the tt 
process is considered uncorrelated and therefore independently altered. 

The measurement of the electron Ej has a relative uncertainty of 2.5% in the endcap region, 
while the uncertainties for the barrel region and for muons can be neglected. These variations in 
the electron energy scale are also propagated to the missing transverse energy. This component 
of the #t uncertainty is treated as uncorrelated with the other ^uncertainties. 

Correction factors have been applied to match the trigger-selection and lepton-selection effi- 
ciencies in simulated samples with those in data. These factors are obtained from data by 
using decays of Z bosons to dileptons. In the electron+jets channel, the correction factor is 
0.933 ± 0.032. In the muon+jets channel, it is 0.965 ± 0.004. The uncertainty on the tt cross sec- 
tion measurement due to the uncertainties of the correction factors is evaluated by weighting 
all simulation-based samples according to the ±1(7 uncertainties obtained from these studies. 

We evaluated the systematic uncertainty on the tt cross section measurement induced by the 
imperfect knowledge of the PDF of the colliding protons using the CTEQ6.6 [49] PDF set and 
the LHAPDF [50] package. For this purpose, a reweighting procedure is applied to all gener- 
ated samples, in which each CTEQ6.6 PDF parameter is independently varied by its positive 
and negative uncertainties, with a new weight assigned to each variation. The resulting tem- 
plates are used to estimate the impact of variations in the PDFs on our measurement. 

The default samples used for this analysis were produced without any pileup collisions. These 
samples are insufficient for the simulation of data taken in late 2010, when the instantaneous lu- 
minosity was substantially higher than in early data-taking and roughly four to five additional 
collisions per bunch crossing were expected. In order to estimate the effect of these pileup 
collisions, which are present in data but not in our fit templates, additional samples of tt and 
W+jets events that included the simulation of these extra collisions were produced. Although 
the average number of pileup collisions in these samples is slightly larger than the expected 
number in data, these simulated events can still be used to provide a conservative estimation 
of the impact of pileup collisions on our measurement. 

6.3 Electron+Jets Analysis 

In the electron+jets channel, we model tt, W+jets, Z+jets, and single-top-quark production 
using the simulated event samples described previously. For the QCD multijet background, 
a sideband method based on data is used to model the M3 and Uj distributions, where the 
sidebands are chosen to be depleted in contributions from real W bosons. In the sideband 
selection, events must have an electron that fails at least two of the three quality requirements: 



6.3 Electron+Jets Analysis 
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7 re i < 0.1 (but 7 re i < 0.5 is always required), transverse impact parameter < 0.02 cm, and the 
standard electron identification criteria. As verified with simulated events, the data sample 
extracted in this way has a QCD multijet purity larger than 99%. In addition, the M3 and 
#t shapes derived from this sample are in good agreement with the distributions from the 
simulation. The fraction of events in the three-jet and inclusive-four-jet sample for each process 
are taken from the simulation. Figure [2]shows the distributions of #xand M3 from the simulated 
three-jet and four-or-more jet samples, respectively, for the different processes. 



CD 
> 
CD 



0.3 



0.25- 



_CMS Simulation ats/s = 
e+iets, N =3 

jets 



0.2- 



0.15 



— tf 
Single-Top 

-W->lv 

— Z/y*-»lT 
— QCD/y+jets 




0.05 



CD 



0.25 



0.15 



- CMS Simulation atVs = 7 TeV 

— tt 



~e+iets, N > 4 

jets 



140 160 

£ T [GeV] 




Single-Top 

-W->lv 

— z/y*^i + i' 

— QCD/y+jets 



0.05 



100 



300 400 



500 



600 700 

M3 [GeV] 



Figure 2: Simulated distributions from electron +jets events of (left) for events with three jets 
and (right) M3 for events with four or more jets. The contributions from the different processes 
are shown separately, and are normalized to unity. Error bars are statistical only. 

The modeling of QCD multijet events from data might induce an additional source of system- 
atic uncertainty. This is investigated by separating the sideband region from which the QCD 
templates are derived into two parts. The sideband region is defined, in addition to other cri- 
teria, by 0.1 < J re i < 0.5 for the electron+jets channel. The QCD template is further split into 
two separate samples of equal-width regions in l te \ (0.1 < I re \ < 0.3 and 0.3 < J re i < 0.5), and 
the templates from these two samples are used to estimate this systematic uncertainty. The 
uncertainty in the ratio of the number of events with three jets to that with four or more jets is 
investigated as well. While this ratio for the model predictions is taken from the simulation, the 
observed ratio in the sideband selection is different. Consequently, the two potential sources of 
systematic uncertainties are studied separately via two independent strength parameters. 

The tt cross section is measured, a ccou nting for statistical and systematic uncertainties, using 
the fit method described in Section 



6.1 



The parameter jS^, which is used to compute the tf cross 
section, and the values of fa for the background processes are determined in the fit. The results 
for /3 t t and the signal and background event yields for the inclusive three-jet bin are given in 
Table [2] While the number of fitted tf events N t { is quoted with its combined statistical and 
systematic uncertainty, for the remaining processes only statistical uncertainties are given. A 



list of all systematic uncertainties in this channel is provided in Section 6.5 with the dominant 
systematic uncertainty coming from the lack of knowledge of the jet energy scale. 



In the electron+jets channel, the resulting tt production cross section is: 



(Tjj = 180±H (stat. + syst.) ± 7 (lumi.) pb . (5) 

The fit produces a combined statistical and systematic uncertainty, as given above. A fit using 
only the nominal templates yields a statistical uncertainty of ^22 pb- Assuming uncorrelated, 
Gaussian behaviour of statistical and systematic uncertainties, one can subtract the statistical 
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6 Cross-Section Measurement 



uncertainty in quadrature from the overall uncertainty, resulting in a systematic uncertainty of 



+39 



32 pb. Individual uncertainties are summarized in Section 



6.5 



Table 2: The predicted and fitted values for fi^ and for the numbers of events for the various 
contributions from the inclusive three-jet electron+jets sample. The quoted uncertainties in 
the tt yield account for statistical and systematic uncertainties, while the uncertainties in the 
background event yields are derived from the covariance matrix of the maximum-likelihood 
fit and therefore represent purely statistical uncertainties. 







N tf 


^"single-top 


Nw+jets 


Nz+jets 


Nqcd 


electron+jets (predicted) 


1.00 


325 ± 52 


31 + 2 


468 ± 34 


81 + 6 


367 ± 27 


electron+jets (fitted) 


A - 1 -0.24 


Q71 +94 
0/1 -78 


33 + 9 


669 ± 61 


116 + 36 


422 ± 51 



The measured tt cross section, in combination with the background estimation, can be used to 
compare distributions of #rand M3 found in data with those predicted by Monte Carlo simu- 
lation. Figure [3] shows the distributions of the missing transverse energy and M3 as observed 
in data. For comparison, the templates from simulation are normalized to the fitted fractions. 
The deviation visible in the high-M3 region between simulation and data has been investigated 
using pseudo-experiments including statistical and systematic uncertainties. For 10% of the si- 
multaneous fits to Uj and M3 in these pseudo-experiments, the derived Kolmogorov-Smirnov 
(KS) value is larger than the KS value observed in data. Therefore, the observed deviation in 
the M3 distribution is not outside the range of expected fluctuations. 
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Figure 3: Electron+jets channel: Comparison of the distributions in data and simulation of the 
discriminating variables Uj (left) and M3 (right) for signal and background. The simulation 
has been normalized to the fit results. Only statistical uncertainties are shown. 



6.4 Muon+Jets Analysis 

The same analysis method is used to measure the tt cross section in the muon+jets final state. 
#T and M3 are again used as discriminating variables. Shape comparisons for the different 
physics processes are shown in Fig. |4| In the muon+jets channel, the QCD templates for these 
two distributions are derived from data by selecting events in a sideband region enriched in 
QCD multijet events. The relative isolation is required to be between 0.2 and 0.5 for the side- 
band selection, in contrast to the nominal selection, where the muon must have a relative iso- 
lation smaller than 0.05. The gap between the allowed isolation ranges in the two selections 
reduces the signal events contribution to the sideband. Events containing muons with large 



6.5 Combined Electron+Jets and Muon+Jets Analysis 
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relative-isolation values have different kinematics due to the correlation of the relative isola- 
tion with transverse momentum. We therefore restrict I re \ to be smaller than 0.5. The QCD 
multijet purity as measured from simulation is 98.4% in the three-jet sample and 94.3% for 
events with four or more jets. As in the electron+jets channel, the QCD template is split into 
two separate samples of equal-width regions in I re \, and the templates from these two samples 
are used to estimate the systematic uncertainty on the QCD modeling. Apart from the elec- 
tron energy scale, all other systematic uncertainties described in the previous section are also 
accounted for in the muon+jets analysis. A summary of the individual contributions from the 



various sources of systematic uncertainties is provided in Section 6.5 
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Figure 4: Simulated distributions from muon+jets events of (left) ffjfor events with three jets 
and (right) M3 for events with four or more jets for the different processes. The contributions 
from the different processes are shown separately, and are normalized to unity. Error bars are 
statistical only. 

The results for (5$ and the various background yields from the binned likelihood fit to the 
inclusive three-jet muon+jets sample are given in Table [3] Using the method described in Sec- 
tion 6.1 the fitted value j6 t { and the ±lc statistical + systematic uncertainties corresponding to 
the fitted value tfg = 1.07 are determined. The result of the muon+jets analysis is a measured 
tt production cross section of: 



(Tjj = 168+^ (stat. + syst.) ± 7 (lumi.) pb . (6) 

The statistical uncertainty is +17 pb. With the assumption that the statistical and systematic un- 
certainties are Gaussian and uncorrelated, the systematic uncertainty is calculated to be +3® pb. 
Similar to the measurement in the electron+jets channel, the fitted numbers of W+jets, Z+jets 
events and QCD multijet events are found to exceed the predicted values. The KS p-value of 
this fit result has been determined to be 95%. Figure [5] shows comparisons of the distributions 
of #j and M3 between data and simulation, where the simulation has been scaled to the results 
obtained from the binned likelihood fit. 



6.5 Combined Electron+Jets and Muon+Jets Analysis 



The tt cross section is also determined using the method described in Section 6.1 for the com- 
bined electron+jets and muon+jets channel. Simultaneous fits of the $j and M3 distributions 
are performed in both the electron+jets and muon+jets channels. Six fit parameters are used: 
the fraction of tt events (/3 t t), the fractions of the different background processes (/3 t , j8w, and 
j6z), and two distinct fractions of QCD multijet events (jSqcd,? and Pqcd,^)- The use of two fit 
parameters for the fraction of QCD multijet events is motivated by the fact that the sources 
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Table 3: The predicted and fitted values for /3 t t and for the numbers of events for the various 
contributions from the inclusive three-jet muon+jets sample. The quoted uncertainties in the 
tt yield account for statistical and systematic uncertainties, while the uncertainties in the back- 
ground event yields are derived from the covariance matrix of the maximum-likelihood fit and 
therefore represent purely statistical uncertainties. 







N tt - 


N s ingle-top 


Nw+jets 


Nz+jets 


Nqcd 


muon+jets (predicted) 


1.00 


408 ± 64 


41 ±2 


601 ± 43 


57 ±4 


58±4 


muon+jets(fitted) 


1 07+ ' 26 
A,u/ ~0.24 


407+IO6 


41 ±12 


813 ± 59 


76 ±22 


123 ± 33 
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Figure 5: Muon±jets channel: Comparison of the data distributions of the discriminating vari- 
ables Uj (left) and M3 (right) and the simulation of the different processes. The simulation has 
been normalized to the fit results. 



of such events contributing to this background in the electron ±jets channel are very different 
from those contributing to the muon±jets channel. The cross section was determined with the 
same procedure used for the individual electron and muon channels. Figure [6] shows the Ney- 
man construction with all systematic uncertainties included for the combined measurement. 
The fitted f5 t i parameter and the fitted numbers of events for the various background processes 
are summarized in Table HI 




Figure 6: Neyman construction including all systematic uncertainties for the combined mea- 
surement of the tt production cross section in the electron±jets and muon±jets channels. The 
horizontal line indicates the determined value /3 t t = 1.10 from the binned likelihood fit to ob- 
served pp collision data. 



The fitted value corresponds to a measured tt cross section in the lepton±jets channel of 
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Table 4: Predicted and fitted values for /3 t t and for the numbers of events for the various 
contributions in the inclusive three-jet combined electron+jets and muon+jets sample. The 
quoted uncertainties in the tt yield account for statistical and systematic uncertainties, while 
the uncertainties in the background event yields are derived from the covariance matrix of the 
maximum-likelihood fit and therefore represent purely statistical uncertainties. 





ftt 


N t t 


Ngingie-top 


Nw+jets 


Nz+jets 


Nq C d e+jets 


Nqcd jM+jets 


predicted 


1.00 


733 ±116 


72 ±4 


1069 ± 77 


138 ± 10 


367 ± 27 


58 ±4 


fitted 


1 io +0 - 25 

i - iu -0.21 


8061^ 


76 ±22 


1475 ± 86 


184 ± 51 


440 ±44 


113 ±31 



a tl = 173^32 (stat. + syst.) ± 7 (lumi.) pb . (7) 

The statistical uncertainty is 14 pb. Subtracting this in quadrature from the overall uncertainty 
yields a systematic uncertainty of ^29 P D - The fit in the combined channel yields a KS p-value 
of 68% and agrees well with a simple average of the results in the muon and electron channels, 
while correctly accounting for correlations. 

Table |5]gives an overview of the estimated statistical and systematic uncertainties for this com- 
bined measurement as well as for the two channels separately. The different sources of system- 
atic uncertainties are treated as fully correlated between the two channels, except for flavour- 
specific QCD and lepton uncertainties, which are assumed to be uncorrelated. In order to 
estimate the impact of individual systematic uncertainties, Neyman constructions where only 
the specific source of systematic uncertainty under study is accounted for are used. Each result 
indicates the combined statistical and systematic uncertainties of the contribution under study. 

One can see that the largest contributor to the overall systematic uncertainty is the uncertainty 
in the jet energy scale. Combining both channels significantly reduces the statistical uncer- 
tainty in the measured cross section. However, since both single measurements are already 
dominated by systematic uncertainties, the improvement in the total uncertainty of the com- 
bined measurement is relatively small. 

The combined transverse mass of the charged lepton and the Uj is a kinematic variable that 
lacks the discriminating power of the M3 and %i variables for identifying tf decays. How- 
ever, this variable does provide separation between events containing a decaying W boson and 
non-W-boson decays, and thus serves as an independent check of the kinematics of the simu- 
lated samples used in this analysis. Distributions of the transverse mass in the muon±jets and 
electron±jets channels are shown in Fig.[7]for events with three or more jets. Good agreement 
is found between the data and the sum of the signal and background derived from the simula- 
tion scaled to the fit results. The reduced x 2 value from a fit of the data to the simulation is 1.8 
(0.7) in the electron±jets (muon±jets) channel. 

6.6 Cross-checks 

To test the robustness of the result, the tt cross section is also determined in the muon±jets 
channel using four additional methods. In the first method, we use a procedure based on 
counting the number of events with an isolated muon and four or more jets. This method uses 
an event selection slightly different from that described above. Specifically, the jet pj is required 
to be greater than 25 GeV instead of 30 GeV, and the muon is required to have relative isolation 
less than 0.1, compared to 0.05 in the nominal selection. Also the backgrounds from W/Z±jets 
and QCD multijet events are calculated by using the technique of Berends scaling Ij5lj . In the 
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Table 5: Relative statistical and systematic uncertainties in the estimation of the tt production 
cross section in the electron +jets and muon+jets channels, and their combination, assuming 
j6 t t = 1. The total ("stat.+syst") uncertainty is obtained from a Neyman construction, for 
which all sources of systematic uncertainties are taken into account in the prior predictive en- 
sembles. The estimate of each systematic uncertainty ("syst. only") is calculated by assuming 
uncorrelated, Gaussian behaviour of the statistical and systematic uncertainties and subtract- 
ing the statistical uncertainty in quadrature from the total uncertainty. 





electron+jets channel 


muon+jets channel 


combined result 




stat.+syst. 


syst. 


stat.+syst. 
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uncertainty 


only 


uncertainty 
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uncertainty 
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Stat, uncertainty 
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JES 
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+18.7% 
-15.4% 


+20.3% 
-17.6% 
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Figure 7: Transverse mass (Mr) distributions from data and simulation (scaled to the fit results) 
for (left) electron+jets and (right) muon+jets inclusive three-jet samples. 

second method, we measure the tt cross section using a simultaneous fit to the distributions of 
jet multiplicity (Nj e t s ) arid the muon transverse momentum, pjf. The jet multiplicity has been 
shown (Table [TJ to be a powerful variable for separating top from QCD multijet and W+jets 
events. The variable pj^ is an attractive choice because it is not directly affected by such sys- 
tematics as the JES and JER uncertainties. Furthermore, because the muon in either W or tt 
production comes from a W decay, it receives a significant contribution to its momentum from 
the W rest mass, while muons from QCD multijet events receive no such boost. In the third 
method, the tt cross section is determined from a fit to the muon pseudorapidity distribution in 
order to separate the top signal from the W+jets and QCD multijet backgrounds. This analysis 
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uses the asymmetry between inclusive W + boson and W~ boson production, caused by the 
difference of the quark charges in the initial-state protons, to determine the templates for the 
W+jets background. A fourth method measures the tt cross section from events containing a 
high-px isolated muon and at least three jets. For this analysis, we relax the relative isolation 
requirement to I re \ < 0.1 but introduce a requirement that the #xi n the event is greater than 
20 GeV, in order to keep the amount of QCD multijet background small. A method based on 
Ref . l52l l53l is used to estimate the amount of QCD multijet background separately for events 
with three jets and events with at least four jets. The number of top-pair and W+jets events 
is extracted from a fit to the M3 distribution. All four of these methods give results consistent 
with our previously quoted measurement, but with slightly larger combined statistical and 
systematic uncertainties in each case. 

Complementary to the top-quark-pair production measurements, the cross section for the pro- 
duction of exactly one muon in association with additional hard jets is measured. In all pro- 
cesses considered as signal for this measurement, the muon originates from a W boson. Both 
single top quark decays and decays of top-quark pairs in the lepton+jets channel, including 
decays via tau leptons in the intermediate state, are contributors to this signature. An addi- 
tional component of this signal comes from the production of a W boson with additional jets, 
which is the most prominent background for the analysis of tt "lepton+jets" decays. The same 



event selection as described in Section 6.4 is applied. In addition, all jets in data are corrected 
for pileup, leading to reduced JES and pileup uncertainties. To obtain the cross section, the 
observed number of events in data is corrected for the remaining background processes. These 
include QCD multijet production, the production of a Z boson with additional jets, single-top- 
quark decays, and other tt decays. The number of QCD multijet events is determined from data 
using a template fit to the missing-transverse-energy distribution in each inclusive (or exclu- 
sive) jet-multiplicity bin. The normalization and shape of the other backgrounds is taken from 
the simulation. Figure [8] shows the cross section for the production of a single muon with pj 
> 20 GeV and \t]\ < 2.1 and additional jets as a function of the inclusive and exclusive multi- 
plicity of jets with pj > 30 GeV within \ tj\ < 2.4. The transition from a phase space dominated 
by W+jets events (in the 1-jet and 2-jet bins) towards a region dominated by the production 
of top-quark pairs (in the 4-jet bin) is clearly visible. The comparison of data and simulation 
indicates a good understanding of this transition, while the overall normalization seems to be 
slightly underestimated. This is consistent with the main analysis, which also found a W+jets 
cross section larger than the theoretical prediction. 



7 Discussion and Conclusion 

A measurement of the cross section for top-quark pair production in proton-proton collisions 
at a centre-of-mass energy of 7 TeV has been performed at the LHC with the CMS detector. 
The analysis uses a data sample corresponding to an integrated luminosity of 36 pb -1 and is 
based on the reconstruction of final states containing one isolated, high transverse-momentum 
muon or electron and hadronic jets. The measured cross section for the combined electron +jets 
and muon+jets channels is 173^32 (stat. + syst.) ± 7(lumi.) pb. This measurement agrees 
with but has a larger uncertainty than current theoretical values |lOT[T2l l33l l34j , which agree 
among themselves. For example, the approximate NNLO calculation from Ref. fi"2| yields 
163tl (scale) ± 9 (PDF) pb, while a similar calculation performed using HATHOR (TO! El 
yields 160+^ (scale) ± 9 (PDF) pb. For this calculation, Q 2 = (173 GeV) 2 is chosen for both the 
factorization and renormalization scales and the MSTW2008 NNLO |35j PDF set is used. The 
scale uncertainty is evaluated by independently varying the scales by factors of 4 and 0.25, and 
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Figure 8: Cross section for the production of an isolated muon originating from a W boson 
(including decays via tau leptons in the intermediate state) in association with additional hard 
jets as a function of the (left) inclusive and the (right) exclusive multiplicity Nj e t s of jets with 
px > 30 GeV in the visible range of the detector. The inner error bars on the data points 
correspond to the statistical uncertainties while the full error bars correspond to the statistical 
and systematic uncertainties added in quadrature. The data points are compared with the 
expectation from the event generators used for the simulation. The scaling factors derived in 
the main analysis are not applied. 



the PDF uncertainty is calculated using the 90% confidence level envelope of the PDF pOl [TT| . 
Our cross section measurement also agrees with the earlier CMS measurement in the dilepton 
channel [14] and the ATLAS measurement in the combined dilepton and lepton+jets chan- 
nels |T5| , but has a smaller uncertainty than either of these previous results. Given the agree- 
ment between theory and the experimental measurements in both the dilepton and lepton+jets 
channels, no sign of new physics has emerged in these studies, and the top quark at the LHC 
remains consistent with being a very massive particle whose properties are as predicted by the 
standard model. 
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